Biot-Savart correlations in layered superconductors. 
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We discuss the superconductor to normal phase transition in an infinite-layered, type-II superconductor in 
the limit where the Josephson coupling between layers is negligible. We model each layer as a neutral gas of 
thermally-excited pancake vortices and assume the dominant interlayer coupling is the electromagnetic inter- 
action between the screening currents induced by these vortices. Our main result, obtained by exactly solving 
the leading order renormalization group flow, is that the phase transition in this model is a Kosterlitz-Thouless 
transition despite being a three-dimensional system. While the transition itself is driven by the unbinding 
of two-dimensional pancake vortices, an RG analysis of the low temperature phase and a mean-field theory 
of the high temperature phase reveal that both phases possess three-dimensional correlations. An experimental 
consequence is that the jump in the measured in-plane superfluid stiffness, a universal quantity in 2d Kosterlitz- 
Thouless theory, receives a small non-universal correction (of order 1% in Bi2Sr2CaCu20s+a:). This overall 
picture places some claims expressed in the literature on a more secure analytical footing and resolves some 
conflicting views. 



I. INTRODUCTION 

In this paper, we revisit the problem of the superconductor- 
normal phase transition of a layered, type-II superconduc- 
tor in the limit where Josephson coupling between layers is 
negligible. We model this system as an infinite stack of su- 
perconducting planes, each layer containing a neutral gas of 
thermally-excited, two-dimensional pancake vortices»i We 
take the viewpoint that the dominant mechanism coupling the 
layers is the long-range electromagnetic interaction 2 between 
the screening currents induced by these vortices. We expect 
this model to be relevant to layered superconductors where 
the dominant mechanism by which the superconductivity is 
lost (as the temperature is raised) is the loss of long-range or- 
der in the phase of the order parameter^ Candidate materials 
include the underdoped high-T c cuprates as well as layered 
structures made from conventional type-II superconductors. 

One motivation for considering this model is that investi- 
gations of the superconductor-normal phase transition in the 
cuprates have revealed 3dXY critical exponents, a hallmark 
of Josephson coupling between the planes, in only one case: 
optimally-doped YBa2Cu306+ x , the least anisotropic of 
these materials^ In underdoped B^S^CaQ^Os+z, the most 
anisotropic of these compounds, signatures suggestive of a 
Kosterlitz-Thouless (KT) transition have been seen 5A7 while 
recent measurements 8 on underdoped YBa2Cu306+:r, which 
is substantially more anisotropic than the optimally-doped 
material though less so than underdoped Bi2Sr2CaCu20s+ x , 
show a transition that is neither KT nor 3dXY nor any ob- 
vious interpolation in between. 37 While these observations 
do not conclusively show that Josephson coupling can be 
neglected*^, they do suggest that investigations of different 
mechanisms for coupling the layers might be a fruitful line 
of attack if there is reason to suspect that Josephson coupling 
is very small. The Biot-Savart interaction between screening 
currents in different layers is a long-range, three-dimensional 
coupling which is always present though its influence is usu- 



ally assumed to be small compared to the Josephson term. 

A second motivation is a recent experiment^ on 
La2- x Ba x Cu04 near its stripe-ordered state at x = 1/8 
where a 2d superconducting phase has been observed below 
an apparent Kosterlitz-Thouless transition temperature. One 
proposal 12 suggests that under the right circumstances, the su- 
perconducting state can occur with a finite wave vector, where 
the periodicity is in the same direction as the charge order 
but with double the period. Since the stripes in adjacent lay- 
ers orient with a relative angle of 90°, an orthogonality argu- 
ment implies that the Josephson coupling between first, sec- 
ond, and third neighbor planes is cancelled exactly and further 
estimates 12 imply the residual terms are extremely small. In 
such a scenario, we expect our model to be directly applica- 
ble to the experiments 3 In a similar vein, there are a series 
of somewhat older experiment a 13 ' 14 on various cuprates in a 
dc flux transformer geometry where the observation of effec- 
tively 2d vortices was explained as the cutting of 3d vortex 
lines. Since Josephson coupling is the strongest reason why 
pancake vortices tend to form 3d stacks/lines (as the penalty 
for breaking a line of pancakes bound by the Josephson inter- 
action would be proportional to the system size L) one way to 
interpret the experiment is that the Josephson coupling is ef- 
fectively "cancelled" by the applied field. Once again, in such 
a case we would expect our model to apply. 

Our main result is that the superconductor to normal phase 
transition in a layered superconductor with an infinite num- 
ber of electromagnetically-coupled layers is still a Kosterlitz- 
Thouless transition. The mechanism resembles the single- 
layer problem in that (a) the transition occurs through the 
unbinding of two-dimensional pancake vortices and (b) the 
screening within an individual layer at temperatures above 
Tkt 1S not significantly different from an isolated two- 
dimensional system. However, both the low and high tem- 
perature phases have three-dimensional correlations and the 
jump in the in-plane superfluid stiffness, as inferred from 
a penetration depth measurement, receives a small non- 
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universal correction (of order 1%). Our results are obtained 
via a renormalization group study of the low temperature 
phase and a Debye-Huckel mean field theory of the high tem- 
perature state. 

The electromagnetic coupling may be formulated as an in- 
teraction between the pancake vortices . 1 ' 2 i 15 ' 16 i 17 i 18 ' 19 The ba- 
sic mechanism is that a vortex in one of the layers induces 
screening currents in the same and in other layers, which cause 
Biot-Savart forces on the other vortices. For a single, lit- 
erally two-dimensional layer, this vortex-vortex interaction 
is screened at distances larger than the magnetic penetration 
depth A. For distances much less than A, but greater than 
the coherence length £, the interaction energy of two vor- 
tices of the same sign will be repulsive and scale logarithmi- 
cally with separation. However, if the layer has a small but 
nonzero thickness d, the effective screening length becomes 
A' = 2X 2 /d ^> XJM2. A' often exceeds the sample sizes 
considered in experiment and in such cases, the interaction 
is effectively logarithmic. 

For a layered system, the relevant screening length becomes 
A = 2A|/s, where Ay is the in-plane magnetic penetration 
depth and s the layer spacing. However, for an infinite num- 
ber of layers, the interaction between two vortices of the same 
sign in the same layer is logarithmic at all length scales, not 
just for separations smaller than A. The difference stems from 
the fact that in the infinite layer problem, currents in the other 
layers guide a vortex's magnetic flux radially out to infinity 
within a disk of thickness Ay while in the single layer prob- 
lem, the flux spreads over all spaced For two vortices of the 
same sign in different layers, it turns out the interaction is also 
logarithmic at large distances but attractive Ml Therefore, the 
interlayer coupling favors pancake vortices of the same sign 
aligning into stacks. This is qualitatively what happens with 
Josephson coupling except the attractive force keeping the 
pancakes aligned is now logarithmic and long-ranged, instead 
of linear and short-ranged. 

Because the interlayer interaction is logarithmic, it was 
conjectured^ that the phase transition should be in the 
Kosterlitz-Thouless universality class despite being a three- 
dimensional system. A number of important renormalization 
group studie a 21 ' 22 ' 23,24 , emphasizing the role of pancake vor- 
tices, explored the issue in greater detail. In each case, the 
phase transition was investigated through numerical studies 
of the resulting flow equations, where the interlayer interac- 
tions were treated at varying levels of approximation and de- 
tail (in each case, the Kosterlitz-Thouless equations occurred 
as leading terms when the interlayer interaction was treated 
perturbatively). Refs. pHl22ll23ll provided support for the KT 
scenario while Ref. ||2441 , which was prima facie the most com- 
plete study, reached a very different conclusion: runaway RG 
flows in a very narrow temperature range close to Tkt ap- 
peared to signal a three-dimensional critical regional or per- 
haps a first-order transition. This motivated the present work. 
Our central result, obtained by exactly solving the leading or- 
der RG flows which occur when the interlayer interaction is 
treated non-perturbatively, is that the phase transition is, in- 
deed, in the Kosterlitz-Thouless universality class. 

There are a number of reasons, in addition to three- 



dimensionality, why a KT transition is not a foregone conclu- 
sion for this model. In Ref. 12211 . it was noted that while the in- 
terlayer logarithms individually come with much smaller co- 
efficients than the in-plane logarithm, the infinite set of cou- 
plings obey a "sum rule" (see Eq. ( 12.5b ) which follows from 
flux conservation. Therefore, there is the possibility that the 
collective effect of a large number of layers could influence 
the critical properties. In a similar vein, in Ref. 12511 it was 
noted that configurations involving stacks of vortices where 
one or both ends terminate inside the superconductor, which 
are topologically forbidden in the presence of Josephson cou- 
pling (since vortex lines may then only terminate on the sur- 
face of the material), should be accounted for in its absence. 

In the next Section, we review some basic facts about the 
Biot-Savart interaction. In Section[ni]we present an analytical 
theory of the phase transition using an extension of the 2d mo- 
mentum shell renormalization group 26 . Accounting for con- 
figurations involving stacks of vortices, as discussed above, 
we re-derive the coupled set of flow equations obtained in 
Ref. 1I24I1 and establish this set as an accurate description of the 
low temperature physics. We explicitly solve this set and find 
that the phase transition is, in fact, in the Kosterlitz-Thouless 
universality class. In Sections ITVl and IVl we discuss the low 
and high temperature phases, the latter using a Debye-Huckel 
mean field theory. We conclude with a summary. Technical 
aspects of the calculation are discussed in three appendices. 



II. THE BIOT-SAVART INTERACTION 

In this Section, we give a physical discussion of the Biot- 
Savart interaction and introduce our model Hamiltonian. To 
make our assumptions clear, in Appendix [A] we review how 
this interaction formally arises from a Ginzburg-Landau type 
free energy functional. 

We model our system as an infinite stack of superconduct- 
ing planes, where the stacking is in the z-direction and the 
positions of the planes are given by z n = ns where n is an 
integer and s the interlayer spacing. A two-dimensional "pan- 
cake" vortex 2 ^ of strength mi (mi is an integer) placed at 
the origin of layer n = will induce an azimuthal screen- 
ing current miK c j > (p,ns) in layer n where p is the cylin- 
drical radial coordinate. A vortex of strength ni2 located 
at position (p, ns) will feel a radial Lorentz force given by 
F p {p, ns) = K<t,(p, ns)m 1 m 2 4>o/c where O = §f ■ 

These screening currents were computed in Ref. 1 in the 
limit where s is small compared to the in-plane penetration 
length A|| and for distances p large compared to s and the in- 
plane coherence length £11 : 



K^p,ns = Q) = -^- r 1-— 1 
O7r^Aj|/0 2A|| 



-P/A|| 



(2.1) 



K^(p, ns ^ 0) 



C0Q.S 

167r 2 Aj|p 



(e A « -e A " ) (2.2) 



Eq. d2.ll ) indicates that to leading order, the in-plane screen- 
ing current K^p, 0) ~ -, This implies an in-plane vortex- 
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vortex interaction potential where vortices of the same sign 
repel each other logarithmically with distance, as in the 2dXY 
model.— If, in addition to p ^> s, £y, we also assume that 
p 3> Ay, \ns\, then Eq. (12.21 ) indicates that the out-of-plane 
screening current K,p(p, ns) ~ — -. This implies that the in- 
teraction between two vortices in different planes also varies 
logarithmically with distance but vortices of the same sign 
now attract. Therefore, in this limit, the interaction between 
two vortices of strengths mi and m-2 at positions (xi,nis) 
and (x2, ri2s) is given by: 



V12 ~ -q 2 m 1 m 2 a\ ni ^ n2 \ In 



xi - x 2 



(2.3) 



where 



<*n«{ -*M " .„ ,„ (2.4) 

2A 



s e A V n| ifn^O 



superconducting phase, we expect the density of (thermally- 
excited) pancake vortices in each layer to be small and hence 
the characteristic vortex-vortex separation p to be large. In 
our analysis, we assume that p is large compared to other 
characteristic lengths, such as A||, which is one of the con- 
ditions for the interlayer logarithmic form to be valid. Also, 
in an infinite layer system, the logarithmic approximation will 
break down for interlayer separations \ns\ f» p. However, the 
length scale associated with the convergence of the sum rule 
(Eq. (|2.5l l) is of order Ay so, if p 3> Ay, we expect the net 
contribution of these farthest layers to be small regardless of 
whether the full (Eq. dA13l )) or simplified (Eq. d2.3l )) interac- 
tion is used. Therefore, in this paper, we will approximate the 
Biot-Savart interaction by its long distance form, which per- 
mits us to use Coulomb gas techniques to analyze the partition 
function. However, we will retain some aspects of the short 
distance physics, i. e. on length scales shorter than Ay, by in- 
troducing fugacity variables as discussed in the next Section. 



q = y g^kjii^) 2 ' an d t is a nominal short distance length 
scale (of order Ay). 

Eq. ( 12.3b is a model of the Biot-Savart interaction which 
was considered in Refs.l22l andl24l and is the model consid- 
ered in this paper. A noteworthy feature of this interaction is 
that the coupling constants obey a "sum rule": 



E< 



= 



(2.5) 



III. RENORMALIZATION GROUP ANALYSIS 

The partition function of a layered superconductor, where 
each layer contains a neutral gas of thermally-excited pancake 
vortices is: 

z = E n<? ,, E e *pH 3 E^). c 3 - 1 ) 



{jv M } M 



M 



This is not an accident but follows from flux quantization 22 
and is related to an important feature of the full interaction 
(Eqs. ( 12.1b . (12.21 ). and (lA13b ): the current distribution of an in- 
finite stack of pancake vortices, one in each layer, is exponen- 
tially screened at distances p ^> Ay ,— in analogy with the well- 
known result for a vortex line in a bulk three-dimensional 
superconductor^ In contrast, for a stack of uncoupled two- 
dimensional layers, the current in each layer would decay ac- 
cording to the Pearl criterion^: as ~ - for distances p < A 

2A 2 

and as ~ 4y for p ^> A, where A = d being the layer 
thickness. 

The electromagnetic interaction causes (same charge) pan- 
cake vortices in different layers to preferentially align into 
stacks, which is phenomenologically what happens with 
Josephson coupling except that now the aligning force is log- 
arithmic and long-ranged, as opposed to linear and between 
neighboring layers in the Josephson coupled case. In particu- 
lar, even though a.\ <C ao, which naively suggests that only 
the in-plane interaction is important, the a's decay very slowly 
and the sum rule indicates that the combined effect of many 
layers could possibly lead to three-dimensional effects. We 
will return to this issue in the next Section. 

We conclude this Section by discussing some of the approx- 
imations inherent in Eq. ( 12.3b . In the high temperature super- 

conductors, typical orders of magnitude * are s « 12 A and 



A,, 



1400 A so the smallness of 



10 2 assumed in the 



derivation is met in practice. In and near the low temperature 



where yk : i = cxp(—f3Ek.i), (3 being the inverse temperature 
and Ek.i the energy cost of creating a pancake vortex of type 
k in layer I, the species label k denoting both strength and 
sign. We assume the layers are equivalent which implies the 
fugacities are the same in each layer, i.e. yk,i = yu- Nk,i is 
the number of type k vortices in layer /. The sum on {Nkj} is 
over layer occupations which satisfy charge (vortex) neutrality 
in each layer. The sum on {c} is over spatial configurations 
of vortices consistent with the set {N^.i}^- Vjj is the vortex- 
vortex interaction which includes a hard-core constraint that 
two vortices in the same layer must be separated by a distance 
r of order £y, the in-plane coherence length. 

At this stage, Vij is the exact vortex-vortex interaction. 
To make analytical progress, we separate Vij mt0 two 

parts, Vj + Si^j Vj, corresponding to the contribu- 

tion from pairs of vortices with in-plane separation p < Ay 
and p > Ay respectively. The latter interaction is given by 
Eq. ( 12.31 ) and in a dilute system will apply for most of the 
pairs (with the caveats mentioned the previous Section). We 
approximate the shorter-distance physics in two steps. The 
possibility of having two vortices in the same plane separated 
by a distance less than Ay can be accounted for by suitably 
redefining the fugacity variables.— To approximate the inter- 
action between two closely spaced vortices in different layers, 
we introduce a new set of fugacity variables {w a b-ij}, where 
Wab-ij = cx-p(—f3E a b-ij), E a b-ij being the interaction energy 
of having a pancake vortex of type b in layer j directly above 
a pancake vortex of type a in layer i. For equivalent layers, 
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Wab-.ij will depend only on \i — j\ (and the strengths a and b). 
With these approximations, Eq. ( 13.11 ) becomes: 

Z = 

e n^n^rE-p^E"^ 

{Nk,i;N abiij } k,l ab;ij {c} i^j 

(3.2) 

where N ab - t ij is the number of pairs of vortices where the first 
member is a vortex of type a in layer i which is directly below 
the second member, a vortex of type b in layer j . The sum over 
{c} is now over those spatial configurations of pancakes con- 
sistent with a set of (charge-neutral) layer occupations {N k j} 
and interlayer patterns { N a b-i 3 }. The sum inside the exponen- 
tial is over pairs of vortices that are laterally separated by at 
least Ay and their interaction Vij is given by Eqs. ( 12.3b - d2~4"l i. 

There is another way of viewing Eq. ( 13.21 ), which is more 
in-line with previous treatments of a vector Coulomb gas. Two 
vortices in the same layer are always separated by a distance of 
at least An, which may then be viewed as the effective "size" 
of a vortex. For a system with an infinite number of layers, 
each having a small but nonzero density of vortices, an in- 
finitely long cylinder of radius Ay perpendicular to and pierc- 
ing the layers will "catch" an infinite number of vortices. The 
stack of pancakes caught by the cylinder may be viewed as 
an extended object. By densely packing the system with such 
cylinders, the configurations of the system may be viewed as 
configurations of these extended objects. The extended ob- 
jects can be labelled by a species index n = (. . . , m, 712, • • ■ ) 
where rij is an integer indicating the strength and sign of the 
vortex occupying layer i of the object in question (for a dilute 
system, most of the entries in n will be zero). We can formu- 
late the problem in terms of these extended objects in which 
case the partition function Eq. ( 13.1b becomes: 

z = Eriyn"E cx p(-^E^^)' < 3 - 3 > 

{n} n { c } i^j 

where X){ n } * s a sum over sets °f extended objects consis- 
tent with charge neutrality in each layer and J2{ c } * s ovel c ^ s ~ 
tinct spatial configurations of these objects. The interaction 
between two of these objects, indexed by rij and ir,, is loga- 
rithmic by construction and given by the sum of the pairwise 
interactions between the pancakes comprising each stack: 

V ni , nj =Y J V nik , njl (3.4) 

k,l 

where k, I are layer indices and V nik . njl is given by Eq. ( 12.31 ). 
Since an extended object is composed of an infinite number of 
pancakes, it will require an infinite creation energy E n and the 
corresponding fugacity y n = cxp(—(3E n ) will be formally 
zero4i However, we can write the energy of an extended ob- 
ject as: 

E n = E E n k ,k + E E n k ni;kl (3.5) 
k kjtl 



where the first term is the creation energy of each pancake in 
the stack while the second is the pairwise interaction energy 
between different pancakes in the same stack. Eqs. ( 13 ,4b and 
( 13.5b show that Eqs. (13.2b and ( 13.31 ) are equivalent ways of 
expressing the partition function. 

Our goal is to determine the phase diagram of the model 
described by Eq. (13.21 ) (or Eq. ( 13.3b ) in the dilute limit where 
each plane contains a small, but nonzero, density of vortices. 
The advantage of separating the vortex-vortex interaction into 
a long-distance logarithmic term, accounting for the short- 
distance physics through generalized fugacity variables, is 
that it permits the use of renormalization group (RG) tech- 
niques developed for studying the two-dimensional Coulomb 
gasi 26 ' 28 ' 30 The procedure requires us to consider a generalized 
version of Eq. ( 13.2b : 

Z = 

El TT N k ,l TT N ab:tj TT Nabe-Mk 
UI»*,I Ll W ab;ij 11 W abc; l3 k 
{ Nk t i ; N a b ; ij ... } k,l Q>b;ij abc;ijk 

n ^^'•••)E ex p(-^E > ^) ) (3.6) 

abcd;ijkl {c} i^j 

In this equation, w abc . ljk = exp(-/3E abc . ijk ) where E abc . ijk 
is the three-body interaction energy of an aligned triplet of 
vortices where the first member is a vortex of type a in layer 
i, directly below the second member, a vortex of type b in 
layer j, which is directly below the third member, a vortex of 
type c in layer k. For equivalent layers, E a b c -ijk will depend 
only on \i — j\ and \j — k\ (and the strengths a, b, and c). 
Similarly, we include fugacity variables corresponding to four 
(and higher) body interactions. These terms may be visualized 
in the extended object picture where Eq. ( 13.51 ) generalizes to: 

E n = ^E rlk , k +^E nkrll]k i + ^ E n k nin m ;klm + ■ ■ ■ 
k k^l k<Cl<m 

(3.7) 

In Appendices|B]and|Cl we present a detailed renormalization 
group treatment of this model. The analysis may be viewed as 
an iterative coarse-graining procedure connecting our model 
with a series of other models with the same critical proper- 
ties. A physical picture of this procedure is illustrated and 
discussed in Fig.Q] 

The analysis of Appendix|C]yields an infinite set of coupled 
flow equations for the extended object fugacities {y n } and the 
analogous equations for the ({y k j}, {i% ; fcz}, ■ • • ) variables. 
The system has a fixed point when the {y n } variables are 
identically zero. A linearized theory about this fixed point 
(Eq. dC21l )) suggests that at sufficiently low temperatures the 
fixed point is stable while at higher temperatures these vari- 
ables become RG relevant. Beginning in the low temperature 
phase and raising the temperature, the phase transition is indi- 
cated by one of the {y n } becoming marginal. Which fugacity 
is the first to "unbind" depends on the initial conditions of the 
RG flow. 

If the starting model is given by Eq. d2.4| ), then the first ex- 
tended objects to become marginal are those where n has +1 
in one of its entries and zeroes everywhere else. By Eq. ( 13.71 ). 
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(a) 



(b) 



(c) 



(d) 




FIG. 1: Pictorial summary of the renormalization group procedure 
discussed in Appendices[B]and[C] If an object and its anti-object are 
very closely spaced, then they will cancel upon coarse graining, but 
in the process the coupling constants will be renormalized. (a) shows 
the simplest such process involving single pancakes which may be 
viewed as integrating out the smallest loops of magnetic flux, (b) is 
an example of the more complicated object/anti-object cancellations 
which correspond to integrating out larger loops. Nearby objects that 
are not antiparticles will fuse together, (c) shows two single strength 
pancakes fusing into a double strength pancake, (d) shows a height 
two stack and a single pancake fusing into a height three stack. 



the fugacity of this object is precisely the fugacity of a single 
strength pancake y (= yi i but since we have assumed a trans- 
lationally invariant system, the layer index is not needed). The 
flow equation for this quantity, in the approximation where we 
keep only the leading y dependence, is: 



de 



y 2 {A-(3q 2 a ). 



(3.8) 



In the same limit, the flow equations for the coupling constants 
are: 

d(/3q 2 a n ) 



de 



t)(/3<7 a r , 



(3.9) 



where we have also taken the "distance to marginality", (4 — 
f3q 2 ao), as a small parameter. We will discuss the flow equa- 
tions for the other variables further below. 

Eqs. ( 13.8b and ( 13.9b are precisely the flow equations ob- 
tained in Ref.|24l using a formulation that involved only single 
pancakes, without accounting for extended objects. Our anal- 
ysis implies that the more complicated objects are irrelevant 
at temperatures below the transition so in this limit, Eqs. ( 13.8b 
and ( 13.9b may have greater validity than initially suspected. 
As discussed in the Appendix|C] this irrelevance does not sim- 
ply follow from the linearized theory but involves using the 
sum rule (Eq. ( 12.51 )) to place a bound on higher order terms - 
a technical issue that does not arise in the single-layer prob- 
lem. 



We continue our analysis by recognizing that the right side 
of Eq. ( 13.91 ) is a convolution of the couplings. Taking the 
Fourier transform, we obtain: 



d(l/((3q 2 a(k))) 2 
de 



(3.10) 



where a{k) = J2 n ot n e~ and we used the fact that a m = 
a_ m . Because the right side is independent of k, we may 
formally integrate this equation to obtain: 



[3q 2 a(k, e) 



f3q 2 a{k,0) 



1 + f3q 2 a(k,0)C(e) 



(3.11) 



where C(e) = f Q wy 2 is an integration constant that obeys the 
flow equation: 



dC 
~de 



= Try 



(3.12) 



with the initial condition C(0) = 0. Observe that a(0, e) = 
Yin °- n so ^1- A3. lib implies that the sum rule is preserved by 
the flow. 

Using Eqs. ( 13.8b and ( 13.12b . we obtain a differential equa- 
tion for the flow trajectories in the (C, y)-plane, which can be 
integrated: 



2 2 i 

y =y + 



i 



4C 



71 dk 
.„ 27r 



ln(l+ f3q 2 \a(k,0)C{l)\ 



(3.13) 



Representative trajectories are sketched in Fig. [2] for a given 
value of yo as a function of temperature. As C increases 
from monotonically during the flow (Eq. (13.12b ). the system 
moves along these curves from left to right. At high temper- 
atures, the fugacity increases monotonically during the flow. 
As the temperature is lowered, the fugacity initially decreases 
before increasing. At a critical temperature, the curve will in- 
tersect the y = axis at one point and at lower temperatures, 
the curves cross the axis. However, y = is a fixed point of 
Eqs. (13.8b and ( 13.9b , which means the system will "stop" once 
y = is reached. These curves demonstrate the existence of 
a low temperature phase, where the fugacity renormalizes to 
zero, separated from a high temperature region by a transition 
corresponding to the unbinding of 2d pancake vortices: this, 
by definition, is a Kosterlitz-Thouless transition^ 

The most distinguishing characteristic of the single- 
layer KT transition is the universal jump in the superfluid 
stiffness*^ To see what happens in the layered case, we need 
to determine the range of initial conditions for which the tra- 
jectory defined by Eq. (13.13b passes through y = at some 
C > 0. That is, we need to solve: 



*"Z/o 



f dk 

/ — ln(l + (3q 2 \a(k, 0)|C) - AC (3.14) 

J 27T 



If we begin close to the critical point, and assume the position 
of this point will only be changed by a small amount relative to 
the single-layer case, we can expand this expression treating 



yo, C, and ((3q 2 ao — 4) as small parameters. Then, to leading 
order, we obtain: 



(/3g 2 a (0)-4)C-i^a 2 m (0)C 2 (3.15) 



where we have inverted the Fourier transform. In order to have 
a y = solution where C > 0, the following criterion must 
be satisfied: 



f3q 2 a Q (0)-4> 



,(0) (3.16) 



y 



0.15 




c 



Therefore, at temperatures T > q 2 o;o(0)/4, there is no so- 
lution and y will diverge as the flow coordinate e — > oo. 
Eq. ( 13.111 ) indicates that the couplings {a n } will go to zero 
in this same limit. The critical temperature at which a y = 
fixed point exists is T = g 2 ao(0)/4. This solution represents 
a critical surface (y = 0, {a„ = a„(0)}) where the only con- 
straints on the couplings {a„(0)} are the sum rule (or, more 
precisely, that the matrix aiij = ot\i-j\ is positive definite — 
see the discussion in Appendix iBt and that the values are such 
that the single strength pancake is the first fluctuation to be- 
come marginal. 

To relate this to the universal jump, we need to relate the 
quantity an(0), given in Eq. (12.4t . to the in-plane superfluid 
stiffness measured in an experiment. The in-plane superfluid 
stiffness is defined in terms of the in-plane magnetic pen- 
etration depth Am, which is a directly measurable quantity: 
p s = 0qs/(167t 3 A 2 ) = q 2 /2ir (see the discussion between 
Eqs. (lAlb and (IA2I) ). As the critical temperature is crossed, 
the quantity j3q 2 ao(oo) jumps downward from 4 to zero. In 
terms of p s : 



p s a (oo) 
T 



T+ 



2 

7T 



1 2_ 

1 2A,| 



(3.17) 



(3.18) 



Therefore, while the jump described by Eq. ( 13.17b is a univer- 
sal quantity, the jump in the superfluid stiffness (Eq. ( 13.181 )), 
which is the quantity that is directly measured, receives 
a non-universal correction on the order of 1 percent for 
Bi 2 Sr 2 CaCu 2 8+;c . 



IV. LOW TEMPERATURE PHASE 

Having established the existence of a critical point and 
phase transition, we now turn to the nature of the low temper- 
ature phase and the way the present case differs from a stack 
of decoupled layers. The most direct difference follows from 
Eq. ( 13.1 II ). As the low temperature phase is characterized by 
a finite value of C(oo), Eq. (13. lit indicates that the interlayer 



FIG. 2: Plot of the flow trajectories of Eq. < |3. 13b for the case where 
the initial fugacity yo = 0.1. The six curves, from top to bottom, are 
for f3q 2 = 1, 2, 3, 4, 5, and 6 respectively. Since the flow stops when 
y = 0, this initial condition implies a phase transition occurring for 
[3q 2 between 5 and 6. 



couplings {a n } will have nonzero fixed point values. There- 
fore, interlayer correlations will be present. 

Another perspective may be gained by considering the flow 
equations for the pair interaction fugacity for single strength 
pancakes in layers 1 and m, which we label W\ m (Eq. dC25l )): 



dwi n 

de 



[2 + /3g 2 a TO _i] Wi m - (iq 2 a m -i 



(4.1) 



where we have dropped terms that flow to zero in the low tem- 
perature phase as e — > oo. Eq. ( 14.11 ) has a fixed point when: 



w\ m {oo) 



(3q 2 |a m _i(oo)| (3q 2 |a m _i(oo)| 



2 + Pq 2 |a m _i(oo) 



(4.2) 



In the decoupled limit, u>i m (oo) would be zero. The y = 
fixed point means that at the longest length scales, the ef- 
fective model is equivalent to one where the vortices are not 
present: i. e. the properties of the superconductor will be gov- 
erned by the spin wave part of the action (Eq. ( 1A1 lb ). How- 
ever, at intermediate length scales, which correspond to the 
points forming the RG trajectory, the effective model will be 
one where the vortices are present, albeit with small fugac- 
ity. wim flowing to zero means that configurations where vor- 
tices lie on top of one another contribute progressively less to 
the partition function, relative to configurations where they do 
not, as the system is viewed from larger length scales. The 
reason is purely entropic4^ Eq. ( 14.21 ) indicates an enhanced 
probability for such configurations in the coupled layer case. 

At one level, this is not so surprising because at the small- 
est length scales, the interlayer couplings (Eq. ( 12.41 )) favor the 
formation of stacks. However, it is interesting to discuss the 
higher-body terms. For example, we may obtain the three- 
body interaction fugacity for single strength pancakes in lay- 
ers a, b, and c, which we label w a b c - To obtain this quantity, it 
is easiest to start with the flow equation for y n = y a t, c where 
n is a vector with +1 in layers a, b, and c, and elsewhere. 
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From AppendixICl it follows, to leading order in y, that: 



dy a bc 
de 



(3q 2 

2 — (3a + a. ab + a ac + a bc ) 



Vabc 



Vabip-ac + a-bc) + y ac (a a b + a bc ) + Vbc{a a b + a ac ) 

(4.3) 



where y^ is the fugacity of the extended object with single 
strength pancakes in layers i and j, the other layers being 
empty, and = QSk-j-i. We can determine the low tem- 
perature fixed point value of y a b c by setting the right side of 
Eq. ( 14.3b to zero. By Eq. (13.71 ), these fugacities may be un- 
derstood as a product of terms associated with the creation 
and interaction of the stacked vortices: y ab = y 2 w a b and 
y a bc = y 3 w a bW ac w bc w a bc = y 3 v a bc- Using Eq. (072), we 
find: 



Vabc(oo) 

_ [(<*«& + aac + Ubc) 2 ~ {a 2 ab + a lc + "be)] 



( 3/3q 2 



a - 2) + S±. (a ah + a ac + a bc ) 
@f [(a o6 + a ac + abc) 2 - (a 2 ab + a 2 ac + a 2 c )] 



(^ao - 2 



(4.4) 



In terms of the variable w a bc, the result is: 

«w(oo) 

(w) ^ aab + aac + abc ^ ~ ( a "b + a ac + a lc)] 



aabUacOtbc (^2~ a Q ~ 2 



(4.5) 



In decoupled limit, v a b c and w a b c flow to zero and infinity re- 
spectively, which again is purely entropic. v a b c being nonzero 
indicates an enhanced probability for these objects in the cou- 
pled layer case. However, these expressions indicate that the 
effective interaction between the vortices within the stack is 
no longer a simple pairwise form. Qualitatively, this is another 
indication of three-dimensional correlations being present in 
the low temperature phase. 



V. HIGH TEMPERATURES 

Fig. |2] shows that the high temperature phase has C — > oo 
as e — > oo. From Eq. ( 13.1 lb , we conclude that the new fixed 
point model is one where all of the couplings have renormal- 
ized to zero. In the same limit, y also diverges, which may 
be interpreted as an unbinding of pancake vortices. There- 
fore, the high temperature phase may be viewed as a stack of 
nearly independent planes, each containing a neutral plasma 
of weakly interacting vortices. However, once y becomes of 
order unity, the renormalization group approach discussed in 
the previous Section no longer applies. The flow equations 



(13.8b and (13 - 9b were obtained by dropping terms involving 
higher powers of y, which is no longer valid when y is order 
unity. Physically, the RG approach of Appendix|C] assumes a 
dilute gas of vortices, which is no longer true at high temper- 
atures. 

Therefore, to gain insight into the nature of the high tem- 
perature phase, a different analytical approach is needed. We 
study this limit using a Debye-Huckel mean field analysis 
where we assume each layer has N positive and N nega- 
tive vortices, which we represent via density functions (x) 
and p^(x), where m is a layer index and x is the in-plane 
coordinate. The total charge density in layer m is given by 
p m (x) = /^(x) — p~(x). The system is modeled by the 
mean-field free energy functional: 



d 2 xd 2 yp m (-x)p n (y)V mn (y. - y) 



d 2 x pt (x) ln(^)) + p~ (x) ln(^) 



iV 



N 



J d 2 xc l (x)p m (x) 



(5.1) 



where the three terms are the vortex-vortex interaction, sys- 
tem entropy, and interaction of the vortices with an exter- 
nal potential <^*'(x). The potential V mn (x — y) is given by 
Eq. ( 12.31 ). The idea is to minimize F with respect to the func- 
tions /9,^j(x) subject to the constraint: J d 2 xp^ l (x) = N. 

The previous paragraph is a natural way to discuss the high 
temperature limit of a layered Coulomb gas interacting via 
Eq. (12. 31 ). 32 However, there is a caveat to note when we relate 
this model to superconductors. The phase transition discussed 
in the previous Section describes a loss of superconductivity, 
which is the low temperature state, due to a loss of long-range 
order in the phase of the order parameter. Therefore, in prin- 
ciple, one may have a regime without superconductivity but 
where the amplitude of the superconducting order parameter 
is nonzero^ It is in such a regime that our model applies since 
it is reasonable to assume the basic degrees of freedom would 
still be vortices. A nonzero order parameter amplitude means 
the effective penetration depth is finite so we expect the inter- 
layer mechanism of Eq. (12.3b . which is ultimately due to the 
screening currents, to still apply. On the other hand, once the 
order parameter amplitude is zero, we can no longer think of 
the basic degrees of freedom as vortices so our model would 
no longer be relevant. In Ref. yfl, it was suggested that for 
layered superconductors with small superconducting carrier 
densities, including the high-X" c materials, the energy scale 
associated with phase decoherence might be appreciably less 
than the energy scale at which the amplitude goes to zero. 
In such a case, there would be a temperature range above T c 
where our model would apply. 

In the absence of an external potential, the minimum 
free energy is obtained when both the positive and negative 
charges are uniformly distributed in each layer, i.e. p„(x) = 
N/A = po where A is the area of a layer. If we perturb the 
system about this limit with a small cjf xt , the corresponding 
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density fluctuation may be calculated as a linear response: 
6p m (x) = - I d 2 x' X mn(x - x')CV) (5-2) 

n •* 

where the (Fourier transform) of the susceptibility 

Xmn(x,x') = -((5 ( o m (x)/^ xt (x'))0=o is given by: 



x(q, k) 



l 



(5.3) 



where ^(q, k) = a(k)^§- is the Fourier transform of Eq. ( 12.3b 



and a{k) the Fourier transform of Eq. ( 12.41 



t;(fc) = q(r 



—ikn 



a + a n e 



2 A i 



r e s A|l 
L/M|| 



Mil - 1 



COS 



fc - 1 



cos k — cosh v- 



(5.4) 



We choose cff xt (x.) to be the potential of a unit test charge at 
the origin of the n = layer: ext (q, k) = |ja(fc) While this 
potential is not small near the origin, our interest is in the long 
distance behavior. The corresponding density fluctuation is: 



Sp(q,k) = - X (q,A:)^ xt (q,fc) 



q 2 + ?/ 2 



cos k — 1 



cos k — a„ 



(5.5) 



77 w k(1 — ^-). However, a more striking difference is that if 
we form an infinite stack by placing a test vortex at the origin 
of every layer, then the sum rule indicates that the stack will 
be completely screened while in the case of decoupled lay- 
ers, a each test charge will influence a density fluctuation in 
its own layer of the 2d Yukawa form. In this sense, the high 
temperature phase of the infinite layer model does not corre- 
spond to a complete layer decoupling but retains some of its 
three-dimensional features. 



VI. CONCLUSION 

In conclusion, we have shown that the superconductor-to- 
normal phase transition in an infinite-layered, type-II su- 
perconductor, in the absence of Josephson coupling but in 
the presence of electromagnetic coupling, is a Kosterlitz- 
Thouless transition. The jump in the in-plane superfluid 
stiffness, which is a universal quantity in the single layer 
problem, acquires a small non-universal correction. We find 
that the phase transition is driven by the unbinding of two- 
dimensional pancake vortices but both the low and high tem- 
perature phases show three-dimensional characteristics. 

A natural topic for future work is to find more connections 
with experiment, including ways to distinguish the electro- 
magnetically coupled problem from the single layer case. A 
more pressing issue would be to explore how these conclu- 
sions are affected by having a small but nonzero Josephson 
coupling and/or other mechanisms of coupling the layers. 



where r\ 



2 47rpo s 

T 2A,i 



- A 



+ 1 



and a q = 



q cosh 



Tak- 



ing the inverse Fourier transform in the layering direction: 



<W(q) = 



>m,0 



'a g — l 



al-l 



(5.6) 



The inverse transform in the q direction is not a simple ex- 
pression. However, if we expand in the small parameter j-, 

and also assume that r) is not too large, then to leading order, 
we obtain: 
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APPENDIX A: RELATION OF GINZBURG-LANDAU AND 
COULOMB GAS MODELS 



-t;|x| 



-i7|x| 



>r 



(5.7) 



At large distances, the in- (out-of) plane fluctuation is posi- 
tive (negative) which is expected by neutrality since a positive 
charge at the origin will attract negative (positive) charges to- 
wards the origin in the same (different) layers. The in-plane 
density fluctuation has essentially the same Yukawa form as 
the single layer problem. The screening length is slightly 

renormalized from the in-plane value of k — J ^p- to 



In this Section, we review how the Coulomb gas description 
of a layered superconductor arises from a phenomenological 
free energy functional of the Ginzburg-Landau type. Variants 
of this derivation may be found in a number of references, in- 
cluding the original paper of Efetov^. Our presentation closely 
follows Ref.El 

Our starting point is the Lawrence-Doniach model of lay- 
ered superconductors^ where the system is modeled as a dis- 
crete set of superconducting layers stacked in the z direc- 
tion. The layers are assumed to have the same thickness d 
and are uniformly spaced with interlayer separation s. With 
each layer n, we associate a superconducting order parame- 
ter ^ n = \^ n \e ie " which we assume does not vary in the z 
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direction within a layer. The Lawrence-Doniach free energy functional is then given by: 



J 



■F(*«,A) = dJ2 / rf'V<5(z-z„) 



a|#, 



l 2 + S*„| 4 



^ It / - e * 



2m:s 2 



2toJ 



(n+l)s 



-ihV\\ An 



8- 



d 3 r(y x A) 2 , 



(Al) 



where the subscripts || and z refer to the two in-plane and 
one out-of-plane coordinates respectively and e* = 2e and 
Tom t = 2m || . 2 are respectively the charge and effective 
masses of a Cooper pair. The difference between Eq. ( IAU 
and the usual Ginzburg-Landau functional is that the order 
parameter fields are only defined within the layers so that the 
kinetic energy associated with the z direction is discretized. 
Note, however, that the magnetic field is defined everywhere 
in space. This approach differs from anisotropic Ginzburg- 
Landau theory, where the order parameter is defined every- 
where and the z direction still has a continuum description. 
We expect Eq. (IA1I ) to accurately describe highly anisotropic 
superconductors, such as the high-T c compounds, but we are 
not aware of a precise way in which to derive Eq. (IA1 b as the 
limit of an anisotropic Ginzburg-Landau model. 

Next, we assume the amplitude of the order parameter is 
constant in each layer, i.e. I^nl 2 = n* = where n s is 
the number of superconducting electrons per unit volume in a 
layer. This assumption clearly breaks down within the core of 



a vortex, which is a region with a radius of order £|| = 2 m*|a| ' 

the in-plane coherence length. For the type II superconductors 
of interest in the present work, £y is small compared to the in- 
plane magnetic penetration depth Ay = ( j^-jpi Y^ 2 , which 

is the other in-plane length scale of interest; here (n s ) = n s - 
is the average number density of superconducting electrons 
over the whole sample volume. 1 Therefore, if the smallest 
length scale we are interested in is of order Am and if we 
further assume that the concentration of vortices is dilute, it 
seems reasonable to neglect amplitude fluctuations. 

The term in Eq. dAll i involving m z may then be viewed as 
a Josephson coupling between the phase variables in adjacent 
layers. For the highly anisotropic materials which motivate 
the present work, m z 3> Tom so we expect the Josephson cou- 
pling to be very small. In this paper, we assume the Josephson 
coupling is identically zero and hence ignore this term. 

With these simplifications, Eq. ( lAlb leads to the following 
effective action: 



J 



ns)(X7\\6 n 



2tt 



(A2) 



where 6 n is the phase of the order parameter in layer n; p s = 

- '^Z d is the 2D superfluid stiffness of a layer; and 4>o = ■ 

The interaction between layers is implicit in the second term. 

The next step is to determine the A which minimizes 
the functional ( IA2I) . Once this is obtained, we can rewrite 
Eq. ( IA2b solely in terms of the order parameter. Taking the 
functional derivatives and imposing Coulomb gauge (V • A = 
0) gives the following equations: 



V 2 A Z 



V An 







(A3) 



where i = 4tt Ps (^-) 2 



A 



ns)(An 



00 

2- 



Vn0 n 



(A4) 

Taking the Fourier transform 



of Eq. ( IA4b gives: 

An(q, k) 



A( g 2 + fc 2 ) 



(A5) 



where q and k are the momenta conjugate to the in-plane (xy) 
and out-of-plane (z) coordinates and: 



<*n(q,fc) = X/ 



- ■/ k i 



1 » / x 1 » / . 27TTO. 

A||(q, ns) = - }_^A||(q,fc+ — — ) 

m 

(A6) 
(A7) 



where Am (q, ns) means that the Fourier transform is only in 
the in-plane direction. We can write the analog of Eq. JA5b 
with k replaced by k + (27tto)/s, where to is an integer. 
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Notice that an (q, k + (27rm) / s) = an (q, fc) and <£m (q, k + 
(2nm)/s) = (f»(q,k). Therefore, summing both sides over 
m, and using Eq. ( |A6t , we obtain the minimizing A: 



A||(qfc) 
A z (q, k) 



n (q, ft) 



(g 2 + fc 2 )(A + L(q, fc)) 

q- yn(q, fc) 

% 2 + fc 2 )(A + L(q, fc)) 



(A8) 



(A9) 



where the last expression follows from the gauge constraint 



and 



L(q,k) 



E 



2irm \ 2 



1 



sinh qs 



2q cos ks — cosh qs 



(A10) 



where the last identity is obtained by (standard) complex anal- 
ysis methods. Substituting Eqs. ( IA8b and JA9b into Eq. ( IA2I ) 
will give an effective action in terms of the order parameter 
itself. After some tedious but straightforward algebra, one ob- 
tains: 



T 



32tt 3 A ^ J (2?r 



d 2 g_ 
) 2 



|q x V||6> K (q) 



<Sm,n|q- V|i0„(q)| : 



rffc [q X V||g m (q)][q X V||g ra (q)]e' fc ( m '") s 
2tt g 2 ( g 2 + /fc 2 )(A + L(q,fc)) 



dfc [q ■ V||fl m (q)][q ■ V||0 n (q)]e ifc ( m - n > a 
2tt q 2 {q 2 + k 2 ){A + L(q,k)) 

I 



1 - 



q 2 A 



k 2 {A + L{q,k)), 



(All) 



The two sets of large "()" brackets indicate that the free en- 
ergy functional decouples into two parts which depend respec- 
tively on the divergence-free and irrotational parts of the fields 
{V||#„}. These correspond respectively to the "vortex" and 
"spin-wave" excitations of the order parameter field, which 
appear in the conventional analysis of the 2dXY models We 
make the usual assumption that the phase diagram is deter- 
mined by the vortex part so will ignore the spin waves from 
now on. 

We make the identification 2irp v (q,n) = q x V||0 n (q) 
where p v (x, n) is the vortex number density per unit area in 
layer n42. After performing the k integration^, the vortex part 
of the free energy is given by: 

Fv = J d 2 xd 2 yp v (x 7 m)p v (y 7 n)V mn {x - y) 



where: 



where the vortex- vortex interaction is given by: 

v (x) - *3- f d 2 q e*- 

and 

W mn (q) 



4ttA 2 J (2tt) 2 q 2 



(A12) 

W mn (q)} (A13) 



ssmh(qs) (G q - (G 2 - l) 1 / 2 )!™-"! 



where G a 



2A 2 g 
cosh(gs) 



(Gg 



(A14) 



- 1)V2 

While this interaction 



s sinh(qs) 

in — Luomuoi i T-r^ . 

1 ^ ' 2A|<? 

has a complicated appearance, our interest is primarily in its 
long distance (relative to Am) behavior, which was obtained in 
Ref.lH 



WW) 



in 2 X 2 



ai m -«i In \x\ 



(A15) 



V0 



-|»l 



2A,i 



(A16) 



and the "~" sign indicates that this result is to leading or- 
der in j-, which is assumed to be a small parameter. For 

a discrete set of vortices, p v (x, n) = m^(5(x — Xj) where 
rrii = ±1, ±2, . . . and Xi are the strength and position respec- 
tively of the ith vortex in layer n. If, in addition to the interac- 
tion energy in Eq. dA15l ), we also assume that each vortex has 
a self-energy associated with its core, which we may represent 
through a fugacity, we arrive at the partition function given by 
Eq. ATp . 



APPENDIX B: RELATION OF COULOMB GAS AND SINE 
GORDON MODELS 

The generalized Coulomb gas model discussed in the text 
may be formulated as a sine Gordon field theory, as in the 
two-dimensional case^i: 

S[{z n }, 4>] = - J2 9 -f f d 2 xV^ ■ V^- 

i,3 ^ 

+ / ^cos[n-0(x)]. (Bl) 

In this Appendix, we establish the equivalence of this expres- 
sion with Eqs. ( 13.3b and d3.6l ). In AppendixICl we analyze the 
phase diagram of this action using the renormalization group. 
Here i, j are layer indices and r is a short distance cutoff of or- 
der the in-plane penetration depth An . The matrix g is defined 
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so its inverse g _1 is related to the coupling matrix for the lay- 
ered Coulomb gas, given in Eq. ( 12.4b : g^ 1 /2tt = (3q 2 a\i_j\. 
The factor z n = 2y n , where y n is the fugacity of an extended 
object indexed by the occupation vector n as discussed in the 
text. The vector <fi = (. . . , fa, fa, . . . ) where fa is a sine- 
Gordon field corresponding to layer i. 

We begin by writing Eq. (IB II ) as a partition function Z = 
J D<fie s ^ Zn }'^ and expanding the cosine terms: 



Z = Z °(I1[ 1 + ^ / ^cos(n-0(x)) 

d 2 xd 2 y cos(n- </>(x)) cos(n-0(y))- 



2!r 



o 

(B2) 



Here Z = J Dd>e s °W where S [4>] = -'Eij^ffVfa ■ 
\7<fij and (. . . }o denotes an average over the full Gaussian 
measure. The Coulomb gas partition function follows from 
calculating these averages. We begin by writing the cosine 
terms as: 



(B3) 



n(x)— ±n 



The averages involve calculating expressions like 



(exp [iJ2a=i n a( x a) • <K x a)])o- Performing the Gaus- 
sian integral: 

( e [iEli%(^)^(x.)]) 

= ex p(-^E iVfcAri ^ fc ( )^( )) 



X! n ak{*a)nbl (x b )((fa(x a )(pi (Xfe))- (fa(0)fa (0))) 



abkl 



(B4) 



where fc, Z are layer indices; a, 6 are indices denoting the TV 
objects {n a } entering the average; TVfc = ^ a "fca is the total 
vorticity in layer k due to these N objects; and the two-point 
function is: 



(fa(x)fa{y)) = g kl l 



d 2 q e 



iq-(x-y) 



(2tt) 5 



(B5) 



Notice that (fa(0)fa(0)) = In ± 
sum in the exponential is: 



Hence the first 



- \ J2 W (fa (0)& (0)> = -\ In -N T g- x N + . . . 
2 * — • 2 r 

(B6) 

where ". . . " are terms subleading in L. The form of the 
interaction (Eq. ( 12.41 )) and the sum rule (Eq. ( 12.51 )) ensure 
that the coupling matrix g -1 is positive definite^ Therefore, 
N T g _1 N > 0, implying the average in Eq. ( IB4b will be zero 
as L — > oo unless N = 0. Therefore, the only configura- 
tions of objects (n a , • ■ • , njv) with nonzero expectation value 
are those that satisfy vortex neutrality in each layer, a stronger 



condition than overall vortex neutrality. If this constraint is 
met, then Eq. dB4l) becomes: 

exp[-- ^2 n ak{^a)n b iMV k i(x a - x b )] (B7) 

a, b, k,l 



where k and I indicate the planes where pancakes with 
strengths n a k and nu and in-plane coordinates x a and re- 
side and: 



Vfci(x a - x b ) 



<<MXa)<Mx b )) - (fa(0)fa{0)) 

|x„ - xj 



9u 



2tt 



In - 



(B8) 



where are terms that are subleading for large |xj — xj |. 

The term we are considering will also have a prefactor. Sup- 
pose the N objects entering the average in Eq. (lB4b include 
N na objects of type n a for a = 1, . . . , TV. There will then be a 

factor of n^Li i z na a I (^n Q ) !] from the Taylor expansion and 
a factor of \j1 N from writing the TV cosines as exponentials. 
Eq. ( IB4I ) is the integrand of a 27V dimensional spatial integral, 
which may be viewed as summing over different spatial con- 
figurations of these N objects. If we choose to write this as 
a sum over indistinguishable configurations, then there will 
also be a factor n_=i (\ ) ' f° r the number of configurations 
(contained in the integral) identical to having the objects {n a } 
at spatial positions {x a }. Combining everything (and drop- 
ping the unimportant factor Zq): we may rewrite Eq. ( IB1I ): 



2 = E 11^/2)^ £ exp(- £ V nita . ) 



(B9) 



M 



>r.i 



where the first sum is over sets of objects {n} satisfying vor- 
tex neutrality in each layer and the second is over indistin- 
guishable spatial configurations {c} of these objects. Com- 
paring Eqs. ( IB7l )- dB9l ) with Eq. (13.31 ), we see that the sine 
Gordon theory is equivalent to our generalized Coulomb gas 
if we make the following identifications: 



Vn = Z n /2 



Pq 2 a\ 



;V2tt 



(B10) 



(Bll) 



as asserted. Using Eq. (13. 7b to expand the fugacities, we see 
that the sine-Gordon model of Eq. (IB 1 b is also identical to the 
partition function in Eq. (13.61 ). 

The final point to note is that Eq. ( IB7b takes the logarith- 
mic form at distances large compared to r but will vanish 
as |x a — Xf,| — > 0. In contrast, the interaction that enters 
Eqs. (13. 3b and (13.61 ) is a hard-core interaction (the particles are 
assumed to be at least distance r laterally apart). As discussed 
m Ref. HI, this is not actually a problem as a sine-Gordon 
model literally equivalent to a hard-core Coulomb gas is pos- 
sible with a slight renormalization of the coupling constants 
cjij that will not affect our calculations (since these corrections 
will manifest as higher order terms in the RG analysis). 
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APPENDIX C: DERIVATION OF FLOW EQUATIONS 



cumulant expansion: 



In this Section, we derive the RG flow equations presented 
in the main text. Our starting point is the action for the lay- 
ered sine Gordon model (Eq. (IB lb ~) which, as shown in Ap- 
pendix|B] is equivalent to the layered Coulomb gas (Eq. ( 13.3b ) 
discussed in the text. We analyze this action using an exten- 
sion of the momentum shell renormalizaton group approach 
discussed in Ref. |2g. The calculation applies in the small fu- 
gacity limit. 

The first step is to write the action (IB II ) as S[{z n }, 4>] = 
So[(f>) + Si[{z n }, <fi], where So is the Gaussian term. The 
fields {(f>i} are written as a sum of fast and slow modes, i.e. 

<t>i = </>i,< + <t>i,> where: 



i,<(x) = [ 



|q|e[<M 



^t,>(x) 



d 2 q 
(2ir) 2 



d 2 q 



>i(q) 



|q|e 



(q) 



(CI) 



(C2) 



Here </>j(q) is the Fourier transform of </>i(x), A ~ - is 
an ultraviolet cutoff, and s = 1 + e is a rescaling param- 
eter. The idea is to integrate over the fast modes to get an 
effective action for 0<(x). The Gaussian term separates to 
give: Z = J X>0<e 5o ^<] / P0 >e s °[*>]e firi K* n >'*<'*>] = 
A J P0 <e s 'o[0<]( e -Si[{^},0<,^>]^ > w here the constant A = 

f P0 > e s °[ < ^ > l will be dropped because it does not affect the 
critical properties of the model. The subscript (>) on the av- 
erage, which we will also drop, denotes that the Gaussian av- 
erage is over the fast modes only. 

We can write the average as ( e Sl ^ Zn ^^) = e s 'i^ 2 '>^* < l 
where the relation between Si and S[ may be expressed as a 



s[ = (Si) + -((sf)-{s 1 y 2 ) + .. 



[ d 2 x(cos(n ■ 0(x))> + 

" m,n 

d 2 xd 2 y (cos(m • 0(x)) cos(n • <p(y))) 
- (cos(m • </>(x)))(cos(n • cj>(y))) 



(C3) 



A convenient way to obtain these averages is with the useful 
fact: 

/p rd 3 a;J(xH(x)\ _ e /d 2 a ; J(x)-0 < ( x ) 



x cxp 



k.l J J 



d 2 xd 2 yJ k (x)(^ > (x) ( j )l , > (y))J l (y) 



(C4) 



where the two-point function is: 



9fc,>(X)<R> 



d 2 a e lq ( x " y ) 



(C5) 

'|q|e[#,A] O") Q 

Using these relations with J(x') = m<5(x' — x), we readily 
obtain: 

(cos[n • </>(x)]} = s-^E^fffe, 1 ^"! cos [ n . ^(x)] (C6) 
and with J(x') = in5(x.' — x) ± imS(x! — y): 



(cos[m • </>(x)] cos[n • 4>(y)]) — (cos[m • 0(x)])(cos[n • 4>(y)]) 

= _ s -£i: w 9M{^rn l+ n kni )^ e -i £ w (m k n,+m, rek )(*k,> M) _ ]_] cos [ m . < ( x ) + n . 0< (y)] 

+ [ e * E w (™*«i+™i**W*.>(*)*i.>(y)) _ 1] cos [ m . < ( x ) _ n . 0<( y )]) 
~ \ {j{52 < < mkni + m i n k){<t>k,> (y)}) (cos[m • 0<(x) -n - 0<(y)] - cos[m ■ 0<(x) + n • 0<(y)]) 

kl 

(CI) 

I 



where in the last line, we retained terms to linear order in e. 

As the ultraviolet cutoff A — > oo, the integral dC5b be- 
comes arbitrarily small except when x ~ y. Therefore, 

we make the approximation: (<f>i,> (x)^- > (y)) f» r 2 5(x — 

-i 

y)(0i,>(O)0j,>(O)) w r 2 (5(x - y)^e and replace the co- 
sine operators with the leading term in their operator product 



expansions. 35 For m ^ n, this implies the replacement: 
cos[m • 0<(x) ±n ■ 0<(y)] ~ cos[(m ± n) • </><(x)]. 

(C8) 

Physically, this means that two closely-spaced vortex stacks, 
when viewed from a distance, appear as a single stack com- 
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posed of pancakes that are fusions of those in the two stacks. 
When m = n, one of the terms in Eq. ( IC7b is of this fusion 
type: 

cos[m • (</><(x) + 0<(y))] « cos[2m • 0<(x)], (C9) 
while the other involves a nontrivial operator identification: 

cos[m-(0<(x)-<£<(y))] « -- ^rmmjV^i ■ V<f>j. 

(CIO) 

Physically, this latter term means that two identical, closely- 
spaced vortex stacks of opposite sign, to leading order, screen 
one another at long distances. However, the effect of this 
screening is a renormalization of the interlayer couplings. 

Using these expressions, we may write an equation for 
S[ [{z' n }] to second order in the fugacities: 



S'l[{z'n},4>< 



E ^Jd 2 xcos[n- 0<(x)]. (Cll) 



where: 



7y 



~ 16tt E Z ™ \ E 



(C12) 



A-/ 



and 



87T 



- m k )mig k i 



id 



f 7T- 4+mZra E^ + m *) m l9kl 



kl 



^numg^ . (CI 3) 



4tt 



A7 



The final step is to restore the ultraviolet cutoff of the original 
problem by rescaling the length. The net effect of doing this, 
to leading order in e, is an additional multiplicative factor (1 + 
2e) on the right side of Eq. dC13b . Combining Eq. ( IC11I ) with 
the Gaussian term So, we obtain an action similar to Eq. (IB lb 
but with renormalized parameters. We can use Eqs. dCl ll >- 
( IC13b to write flow equations for these parameters: 



— = | 2 - — V 
de \ 47r 4^ 



n k nig kl 



m \ fci / 

g~ X! Z n+m^m ( E( n fc + m k )mig kl 1 J . (C14) 

m \ kl j 



and: 



de 



7^: E z ™ ( E m k m i9u J ™i TO i- (C15) 

m \ kl ) 



It is more convenient to write the latter in terms of the inverse 
coupling matrix:— 



d 9~i 



de 



16tt 2 2 ™ ( £ 



m k mig kl 




1 m p m q g q ^ 



(C16) 



Eqs. dC14l )-( IC16l i are natural generalizations of the usual 
Coulomb gas flow equations^ However, as mentioned in the 
text, the fugacity z n of an extended object n can be expressed 
in terms of the creation and interaction energies of the pan- 
cakes forming the stack. This latter formulation is convenient 
because in this picture, the fundamental degrees of freedom 
are individual pancakes, which have finite energy, as opposed 
to extended objects, which have (formally) infinite energy. 
The content of Eqs. (IC 1 4b — (IC 1 6b can be recast in this lan- 
guage via the following relation which follows from Eq. ( |3.7b : 



J rii ,nj ;i,3 



i<3 



n 

i<j<k 



(C17) 

where i, j, k are layer indices and the y and w variables were 
defined in the text. This implies: 



dz r 

Zn 



E 



dy ni ;i 



vrii ,rij 



i<j<k 



II! 



(C18) 



The flow equation for the y variables are obtained by consid- 
ering Eq. dC14b for the occupation vector n a: ; which has only 
one nonzero entry: a strength a vortex in layer i. For a transla- 
tionally invariant system, the fugacity variable will not depend 
on our choice of layer i so y a -i = y a and: 



dy a 
Va 
1 



dz r 



o -1 
= 2 4^ 00 



E- 



^2(aS k ,o - m k )mig k ^ . (C19) 



kl 



Using Eq. dC19b , we can obtain the flow equations for the 
two-body fugacities w a b;ij by considering Eq. (IC14I > with oc- 
cupation vector n a b-ij which has only two nonzero entries cor- 
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responding to vortices of strengths a and b in layers i and j. 



dva + dm ]= _ [2 + ^ g 



Va 



Vb 



1 



8tt^ 1 

m 

x (a^mi^ 1 

z ni ;j +m 



+ 



in 



n f>;j 
^n a .i — m ^1 



I 

+m 2 ni);j ._ m 



■Tl oi ,.;,— m 



(C20) 



In this manner, we can systematically obtain the flow equa- 
tions for each of the many-body fugacities, though the formal 
expressions become increasingly complicated. 

From Eqs. ( IC14b -( IC16l ). we see that the system has a fixed 
point when all of the z n 's equal 0. To probe the stability of this 
fixed point, we consider Eq. dC14b . keeping only the linear 
term for the moment. Converting to the variables in the main 
text via Eqs. (IBIOMbTTT ). we have: 



dy n 
de 



= 2- 



E 



n k nia\ k _i\ 



(C21) 



where the term in square brackets is the coefficient of the loga- 
rithmic interaction energy between two stacks of pancake vor- 
tices characterized by the same occupation vector n. From 
Eq. ( |2.4b . it may be shown that this quantity is always non- 
negative and strictly positive for occupation vectors n that are 
"compact", i.e. whose non-zero entries all occur in a region of 
finite extent. For compact vectors, Eq. (IC21I ) indicates the cor- 
responding fugacities will be irrelevant at zero temperature. 
As the temperature is raised, the first fugacities to become 
marginal correspond to the vectors {n^}, which have a unit 
strength pancake in one layer, the other layers being empty. 
Similar considerations apply for non-compact vectors that are 
sparsely filled. The magnitude of the square bracket term can, 
in principle, be lowered by considering non-compact vectors 
that are densely filled. For example, in the case where we have 
a unit strength pancake in every layer of an infinite system, the 
square bracket term will vanish by the sum rule. However, in 
a dilute system, such objects will not be present in the initial 
model. 

Therefore, in the case which concerns us, the stability of the 
fixed point is determined by the {y ni . 4 }, which via Eq. ( IC17b 
is equivalent to the single strength pancake fugacity yi = y 
governed by the equation: 



g = (2-^o) y + (y 3) 



(C22) 



We can similarly approximate the other flow equations by 
keeping only the leading powers in y. Because we are in- 
terested in the critical behavior, we take the "distance to 
marginality" x = ((3q 2 ao — 4) as an additional small param- 
eter. In such an expansion, Eq. dC 16b becomes: 



d([3q 2 a n ) 
de 



-iry 2 J2(fa 2 a n - P )Wq 2 a p ) (C23) 



The simplest composite object is a pair of single strength vor- 
tices, the first in layer 1 and the second directly above in 
layer m. The flow equation for the corresponding fugacity 
2/nunm = Vim simplifies to: 



dyi 



de 



= [2 - /3q 2 (a + a m -i)} yim -f3q 2 a m -iy 2 (C24) 



The only fusion term surviving the expansion is the simplest 
one: single strength pancakes in layers 1 and m combining to 
form the object y\ m . The equation for the two-body fugacity 
Wii-iu = Wi m becomes: 



dwi n 
~de~ 



[2 + /3g 2 a m _i] Wi m - Pq 2 a r , 



(C25) 



In this manner, we can obtain simplified versions of the full 
flow equations appropriate for the physical limit that interests 
us. 

There is an important technical point implicit in this proce- 
dure. In Eq. ( IC22I ). we have asserted that the terms we have 
ignored are 0(y 3 ). From Eq. ( IC 14b . we can see that an ex- 
ample of such a term is g\xy\ m y- This arises from the fusion 
of the extended object composed of two unit strength posi- 
tive pancakes in layers 1 and m with a unit strength negative 
pancake in layer m. Since we expect y\ m ~ y 2 , this term 
will be of order y 3 . However, there are an infinite number 
of such terms which, in total, could potentially overwhelm 
the "leading term" whenever y is nonzero. This differs from 
the usual single plane Coulomb gas where there are only a 
finite number of 0(y 3 ) processes (such as, for example, the 
fusion of a +2 and —1 pancake as well as terms related to 
the cutoff procedure^). To see this is not a problem, note 
that in the low temperature phase, Eq. dC24b implies that 
yim ~ | Q! m|y 2 /2. The infinite sum converges due to the sum 



mle: Em^l VlmV - V Em/l 



aoy . The argument 



is similar for higher order processes. In the absence of such 
a sum rule, the pile-up of an infinite number of higher order 
terms can overwhelm the leading term and hence invalidate 
the analysis 41 
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However, a jump in the superfluid stiffness characteristic of a 
KT transition has been seen in very thin films of underdoped 
YBa 2 Cu 3 6+ x.- 

More precisely, since the rotational symmetry of the plane is bro- 
ken, we would expect the pancakes to be anisotropic in shape, 
since the coherence length will be different whether one is paral- 
lel or perpendicular to the stripes. Also, the interactions them- 
selves might be anisotropic. However, we expect the leading 
long-distance behavior to be at least qualitatively captured by an 
isotropic model. 

A common way of writing Eq. d3.lt is with the sum over {c} re- 
placed by a multidimensional integral over the coordinates of the 
vortice o 24 ' 28 in which case factors l/(Nk,i)\ are needed because 
the vortices are indistinguishable. 

More precisely, the core energy Ek.i will no longer represent a 
"bare" core energy but a renormalized energy that includes the ef- 
fects of screening due to closely spaced vortex-antivortex pairs. 
This screening will also renormalize the coefficients {a n } of the 
logarithmic long-distance interaction from the values in Eq. <2.4t . 
For a sufficiently dilute system of vortices, we expect these renor- 
malizations to be very small. Therefore, we will still use the values 
in Eq. i2Al as initial conditions for the RG analysis of Eq. l !3.2t . 
There is no inconsistency because any physical quantity will in- 
volve an average where a substantial part of Eq. d3.3 b will appear 
in both the numerator and denominator and hence cancel. 
A coarse-graining step may be viewed as taking the effective core 
size of a vortex from r to r + dr . A vortex in the new model repre- 
sents (1 + ^r) 2 configurations of the old model. This is accounted 
for by renormalizing the fugacity. Now consider an extended ob- 
ject yim, with pancakes in layers 1 and m. In the coarse-grained 
model, this object will still represent (1 + ^-) 2 configurations of 
the old model but merely increasing the size of the pancakes will 
give an extra factor of (1 + ^) 2 . Hence, because yi m = y 2 wi m , 
we can account for this overcounting by renormalizing wi m by 
the factor (1 - ^) 2 . 

To see this, note that a unit strength vortex at the origin is repre- 
sented by the phase field V0(x) = j^y . For a contour C encir- 
cling the singularity, § c X79 ■ dl = 2n. In the continuum case, 
this becomes § c V9 ■ dl = 2n J p„(x)dA However, by Green's 
theorem, we know that § C V6 ■ dl — /(Vx V0) ■ d~A. The iden- 
tification 2np v (x.) — (V x V#) ■ z follows. 
To do the integral, it is helpful to first break the range 
[—00,00] into separate integrals over the ranges [— 3n/s, n/s], 
[— n/s, n/s], [vr/s, 37r/s] and so on. In the integral over [(2a — 
l)n/s, (2a + l)7r/s], make the variable change k — > k + 2ira/ s, 
so that all integrals are from [—iv/s,n/s]. Using Eq. dAlOt and 
the fact that Z/(q, k + 2na/s) = -L(q, k), the separate integrals 
can be combined into one which can be performed by contour in- 
tegration over the unit circle. 

To see this, note that one test of whether a symmetric matrix is 
positive definite is if (a) its diagonal entries are positive and (b) 
the diagonal entry is larger than the sum of the absolute values 
of the other entries in its row (or column). A second test is that a 
symmetric matrix is positive definite if and only if all of the lead- 
ing principal minors of the matrix are positive (i.e. which means 
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the following matrices have positive determinants: the upper left 
1-by-l column, the upper left 2-by-2 column, and so on up to the 
matrix itself). The sum rule, and the fact that the magnitudes of 
the couplings monotonically decrease with interlayer distance, en- 
sures, by the first test, that for a finite number of layers, the matrix 
g _1 is positive definite. As these are the leading minors of the in- 
finite matrix, the second test ensures that the positive-definiteness 
holds in the infinite layer case as well. 



To obtain this equation, it is helpful to start from the infinitesimal 
expression for g: g' = g + £7 = g(l + eg _1 7), where the 
matrix 7 is given by Eq. dC12l >, Then, to leading order in e: g~ 1 = 
(1 - eg~ 1 7)g _1 = g _1 - eg _1 7g _1 - 

This pile-up issue is sometimes relevant in theories of sliding Lut- 
tinger liquids. E. Fradkin, private communication. 



